Donor centers and absorption spectra in quantum dots 
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We have studied the electronic properties and optical absorption spectra of three different cases 
of donor centers, D°, D~ and D 2 ~ , which are subjected to a perpendicular magnetic field, using the 
exact diagonalization method. The energies of the lowest lying states are obtained as function of the 
applied magnetic field strength B and the distance £ between the positive ion and the confinement 
xy-plane. Our calculations indicate that the positive ion induces transitions in the ground-state, 
which can be observed clearly in the absorption spectra, but as £ — > the strength of the applied 
magnetic field needed for a transition to occur tends to infinity. 

PACS numbers: 73.21.-b,71.35.-y 



I. INTRODUCTION 

State-of-the-art techniques in the manufacture of semi- 
conductor devices have made it possible to construct 
quasi-zero-dimensional electron systems or quantum dots 
(QDs), also known as artificial atoms. These small solid- 
state systems have similarities with properties of atoms. 
Their charge and energy spectra are quantized, and also 
they have a full energy level structured. Because,-pf this, 
they are sometimes known as "artificial atoms" EJ, or as 
"designer atoms"! 2 ]. Modern techniques like self-organized 
growth or molecular beam epitaxy, allow scientists to 
fabricate devices as small as 10 nm and yet their size, 
shape and other properties can be controlled in the exper- 
iments. Two-dimensional electron systems (2DES) are 
confined laterally in between two semiconductor layers. 
Usually the materials are GaAs and Gai-^Al^As. The 
comparison between quantum dots and real atoms is very 
descriptive, the nuclear attractive potential is replaced by 
an artificially created confinement. 

There is an increasing interest in the electronic struc- 
ture and properties of donor centers in semiconductors 
and quantum dots in magnetic fields. Our work has been 
focused in studying the cases of the neutral donor center 
(D°), negative donor center (D~), which is formed by a 
neutral donor trapping an extra electron!], and the case of 
a neutral donor which traps 2 electrons (D 2 ~). In these 
quasi-two-dimensional systems electrons are confined to 
the xy-plane and a positive ion is located at a distance £ 
from the plane in the z axis. This model was originally 
proposed by Fox an Larsenu. If we restrict to the condi- 
tion £ = we are dealing with strictly two-dimensional 
donor centers. The greatest interest has been paid into 
the case of the negative donor center D~ because it is the 
one of greater physical interest. Riva et ah considered in 
a previous theoretical worku the D~ system for quantum 
well structures. They found singlet-triplet transitions as 
function of the magnetic field and the parameter £. 

The negative hydrogenic donor center D~ , is the solid- 



state analog of the H~ ion, which consists of a positive ion 
with two bounded electrons. This system has attracted 
considerable attention since the early days of quantum 
mechanics. It was as early as 1929, when BetheO pre- 
dicted the stability of the H~ ion. Nowadays it is well 
known that at zero magnetic field H has only one bound 
state, while when magnetic field increases, the ion binds 
additional states. This problem has also attracted at- 
tention in astrophysics, as the negative hydrogen ion has 
been found to be of great importance for the .opacity of 
the atmosphere of the sun and of similar starso. 

The existence of D~ centers in center-doped 
GaAs/Al^Gai-^As multiple quantuiH wells was first re- 
ported by Huant, Nadja and Etiennea in 1990. A system 
in which two electrons confined to the xy-plane are bound 
by a positive ion is called a strictly two-dimensional D~ 
center. Fox and Larseno called it D~ barrier center. It 
is modeled by generalizing the strictly two-dimensional 
model by retaining electron confinement in the xy-plane, 
but moving the positive ion to a distance £ from the 
plane on the z-axis. In recent years, there has been a lot 
of theoretical work on the negative donor barrierEl The 
behavior of the strictly two-dimensional D~ center in the 
middle of narrow q uantum wells is known in the strong 
magnetic field limitBBB. 

In a recent articleEl Brun et al. have made near-field 
spectroscopy measurements of CdTe/ZnTe single quan- 
tum dot at low temperatures. They found in the pho- 
toluminescence spectrum peaks proving the existence of 
neutral and negatively charged excitons (X and X~ re- 
spectively) and also multi-excitons. Charged excitons 
form by hole trapping by residual acceptors in the barrier 
material. They claim that the signature of the double- 
negatively charged exciton A 2 " is present in the pho- 
toluminescence spectra. These latest experiments have 
motivated us in considering the calculation of the far in- 
frared (FIR) absorption spectra for D~ and D 2 ~ centers. 

D~ centers are one of the simplest many-body systems 
which cannot be solved exactly. This system is very in- 
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teresting in the study of electron-electron correlations in 
a two dimensional electron gas (2DEG), as the screening 
of impurities in two dimensions in known to have very 
drastic effects in the state of the system as, for example, 
in the quantum Hall regime. In our model the separation 
from the dot to the positive ion is expressed in units of 
<jg = fi 2 e/m*e 2 (effective Bohr radius). 



II. THE EXACT TREATMENT 

We consider a model of donor centers like in earlier 
works but generalized for a number of electrons n c . We 
work with electron of effective mass m* moving in the 
xy-plane confined to a parabolic potential and subjected 
to a perpendicular magnetic field. In the numerical cal- 
culations that follow, we have used material parameters 
appropriate for GaAs, i.e., e — 13 and m* = 0.067 m c . 
We have included spin in our quantum dot calculations 
but ignored the Zeeman energy. The confinement poten- 
tial strength is chosen to be the one of GaAs dots, i.e., 
huio = 4 meV and the effective Bohr radius is a|j = 102.8 
A. The Hamiltonian for a donor center containing n e elec- 
trons is given by 
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H = H + H D + t; — / -, — 



where Hd is the donor Hamiltonian given by 
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The term Ho stands for the exactly soluble part of the 
Hamiltonian, i.e., its eigenfunctions and eigenvalues can 
be found analytically. This part is written as 



where the vector potential is in the symmetric gauge, 
A(r) = ^B(—y, x, 0), the single-particle energies and 
wave functions can be expressed explicitly as 

E° nl = (2n+|l| + l)ftn-ifcj c I, (4) 
VwM) = exp(-i/#)i? ni (r), (5) 

where R n i(r) is the radial wave function, given by 

R nl {r) = Cexp[-r 2 /(2a 2 ) - il6]r^ {r 2 / a 2 ) , (6) 

in which C is the normalization constant, 
a = y/h/(m*n), fi = ^lj' 2 +uj2/4, and L*(x) 
is the associated Laguerre polynomial. It is conve- 
nient t o define al so the following two parameters, 
b — \/l + 4wg /oj 2 is a dimensionless parameter and 
£o = \JhcjeB is the magnetic length. The radial and 



orbital angular momentum quantum numbers can have 
the following values 



n= 1,2, 



I = 0,±1,±2,. 



(7) 



This single-particle basis is used by the exact diagonal- 
ization method in order to expand the Hamiltonian H. 
As H is independent from spin, we can add to our basis 
set the quantum number s z . Since the only effect of the 
spin degree of freedom is to introduce exchange terms in 
the mutual interaction part of the many particle Hamil- 
tonian (|l|), in what follows we omit its. explicit notation. 

The exact diagonalization methodtJ'Ej consists in ex- 
panding the total Hamiltonian as a series of a given basis 
and extract the lowest eigenvalues (energies) of the ma- 
trix generated. The better the basis describes the Hamil- 
tonian, the faster will be the convergence. A good discus- 
sion of the convergence and comparison of this method 
was made by Pfannkuchc et alt3 The most common ba- 
sis chosen is the one that describes the Hamiltonian at 
zero order, i.e., the basis found from the exactly solvable 
part of the Hamiltonian, which in our case is represented 
by the discrete basis of eigenfunctions V'nZ- 

When generating the Hamiltonian matrix, the part 
corresponding to Ho will lead to the diagonal ele- 
ments (^) , but the part correspondent to the donor center 
is given by the more complicated integral 
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ip n2 i 2 {r)dr. (8) 



The integration over the angular coordinate 9 gives the 
conservation law of the total angular momentum L = 
y^l. Therefore, there will be no mixing of energy levels 
with different total angular momenta and we will be able 
to classify the eigenfunctions by the L quantum number. 

If we want to calculate the matrix elements which are 
different from zero, we must have 1% — l 2 . Under these 
conditions, we can set l\ = I = 1%. We have obtained 
an exact analytic solution for the integral which appears 
in (||), ending up with the analytic solution for the one- 
particle matrix elements 
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where n = Ki + k 2 + Ul + 1, and the coefficients IC? are 
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FIG. 1: Lowest energy levels for a donor barrier D~ with two 
different values for parameter £ as a function of the magnetic 
field. Upper picture corresponds to £ = ag and lower picture 
to £ = |aJa. Some of the levels are labeled in accordance to 
two quantum numbers \L S } . 



FIG. 2: Optical absorption spectra of a donor barrier D~ 
with two different values for parameter £ as a function of the 
magnetic field. The areas of the filled circles are proportional 
to the absorption intensities. 



defined recursively in the following form 
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We must also calculate the two particle matrix ele- 
ments, i.e., the terms corresponding to the electrostatic 
Coulomb interaction between the electrons, for which an 
analytic expression can be found in Ref. |l7]. 

We have also calculated the intensities of allowed opti- 
cal transitions within the electric-dipole approximation. 
If we define the single-particle matrix elements 



d\\> = (A'|re i9 |A) 



2tt5, 



r 2 R y {r)R x (r)dr, (12) 



where A represents the quantum number pair {n, 1} and 



R\ (r) is the radial part of the wave function correspond- 
ing_to the state A. The dipole operators can be written 



x = \ Eaa' Ma'a + d X y] a\,a\, 
Y = k Eaa' - d AA'] a{,a x . 



(13) 



We will apply this equation for the probability of transi- 
tion from the ground state |0) to an excited state |/). 
It can be written as 



A 



</|r|G) 



(f\X\0)\ 2 + \(f\Y\0)f 



(14) 



The areas of the filled circles in Figs. || and [| are propor- 
tional to this quantity. 



III. RESULTS AND DISCUSSION 



We have made simulations for two different values of 



the parameter £, £ 
D~ and D 2 ~ . 



and £ = jag, for both cases, 
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FIG. 3: Lowest energy levels of a donor barrier D 2 ~ with two 
different values for parameter £ as a function of the magnetic 
field. Upper picture corresponds to £ = a% and lower picture 
to £ = ^cib- Some of them are labeled in accordance to two 
quantum numbers \L S ) . 



FIG. 4: Optical absorption spectra of a donor barrier D 2 ~ 
with two different values for parameter £ as a function of the 
magnetic field. The areas of the filled circles are proportional 
to the absorption intensities. 



Figure |l| shows the lowest lying energy levels for a 
donor barrier containing two electrons as a function of 
the magnetic field. Upper and lower plots correspond 
respectively to £ = ctg and £ = \a^. We can see in 
both cases that, at least, a transition in the ground state 
energy level close to 4 and 7 Tesla occurs. This fact is 
also clearly seen in Figure ||, in which the dipole-allowed 
absorption energies are plotted. There we can observe 
that at zero magnetic field it is equally probable for the 
system to jump to the first and second excited states, 
which are degenerate and have L = ±1. At the point 
where the transition occurs, there is a discontinuity in 
the spectrum, which is due to the change of the ground 
state induced by the donor center, i.e., the spin of the 
ground state changes from S = to S = 1. Similar re- 
sults have been obtained in the FIR absorption spectra of 
a confined 2DEG with a Coulomb impurity using Hartree 
and Hartree- Fock methods .Ejcj 

This behavior is also appreciable in the energy spec- 
trum of two and three interacting electrons without any 
donor impurity, but unlike the previous case, the FIR 
absorption lines do not present any discontinuities nor 



splitting. This is due to the fact that two transitions oc- 
cur simultaneously, in the ground or initial state and in 
the excited or final state. The ground state of the D~ and 
D 2 ~ is strongly dependent to the intensity of the raising 
magnetic field, increasing its angular momentum. The 
higher angular momentum means that the average dis- 
tance of electrons from the center of the dot is increased, 
miminizing the Coulomb interaction between them. The 
same effect is also observed in the X~ magnetoexciton.Eil 

In Figures || and [| the energies and the far infrared 
transition spectra are shown for the case D 2 ~ . We ob- 
serve the same effect than in the case of the two elec- 
trons, i.e., there are induced transitions in the ground 
state (|l-§) —* |2 |) — » 3 |)) as B increases, which is 
reflected in the FIR absorption spectra as discontinuities. 
These spin-singlet-spin-triplet transitions have been al- 
ready observed in quantum dots and are knawn as magic 
magnetic number ground state transitions .E3 In this case 
the absorption spectra is richer and presents more strik- 
ing effects, for example, additional absorption lines ap- 
pear comparing to the case of three interacting electrons 
in a quantum dot. It is convenient to mention that the 
optical absorption spectra for the cases of two and three 
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electrons is exactly the same, this is can be seen as a 
result of Kohn theorem. 

Another effect is the splitting of some of the absorp- 
tion lines due to some broken degeneracies in the energy 
levels. We can also see that the smaller is the parameter 
£ the larger is magnetic field needed for the transition to 
occur. Inspecting the behavior of the transition magnetic 
field and the value of £, we can expect that for £ = an 
infinite magnetic field is needed for a ground state tran- 
sition to happen. 

In our model we used an infinite parabolic potential 
in order to confine spatially the electrons, this is not a 
appropriate selection if we want to calculate the bind- 
ing energies of the donor centers D~ and D 2 ~ because 



there is no possibility for electrons to jump to an un- 
bound sta te . We have set a potential depth according to 
referencesc2lcSE3 which allows roughly six bound states 
in the QD. This gives a binding energy for the D~ of 
Eg~ « 3.8 meV and for the D 2 ~ of E^'~ « -8.2 meV 
at zero magnetic field and £ — ajg. These results are in 
accordance with the actual theories and experimental re- 
sults forr-the,!? - which show the stability of the charged 
complexESEa. The three electron system D 2 " , although 
it has localized states, it has a negative binding energy 
and, therefore, it is unstable toward the dissociatiop-into 
a D~ and an free electron in the conduction band.EH 
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